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^ \ Abstract: We propose a new matrix model describing multi-baryon systems. We 

^ \ derive the action from open string theory on the wrapped baryon vertex D-branes 

O \ embedded in the D4-D8 model of large Nc holographic QCD. The positions of k 

baryons are unified into k x k matrices, with spin/isospin of the baryons encoded in 
a set of fc-vectors. Holographic baryons are known to be very small in the large 't 
' Hooft coupling limit, and our model offers a better systematic approach to dynamics 

^ \ of such baryons at short distances. We compute energetics and spectra {k = 1), 

and also short-distance nuclear force {k = 2). In particular, we obtain a new size 
of the holographic baryon and find a precise form of the repulsive core of nucleons. 
This matrix model complements the instanton soliton picture of holographic baryons, 
whose small size turned out to be well below the natural length scale of the approx- 
imation involved there. Our results show that, nevertheless, the basic properties of 
holographic baryons obtained there are robust under stringy corrections within a few 
percents. 
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1. Introduction 

Nuclear physics is one of the oldest branches of high energy physics, yet remains one 
of more difficult. Despite the fact that we know the underlying fundamental theory, 
i.e. QCD, we are still unable to predict, reliably and analytically, behavior of nuclei 
or even a single proton. The problem is of course that one must understand the 
strong-coupling regime of QCD, which by and large remains inaccessible except by 
large-scale lattice simulations. Traditionally, this set nuclear physics apart from the 
rest of high energy physics in many aspects. However, recent developments in the 
so-called gauge/gravity duality began to solve certain strongly coupled field theories, 
possibly including QCD or its close relatives, allowing the two communities merging 
with each other. In this paper, by using the gauge/gravity duaUty, we propose a 
new matrix model for the dynamics of multi-baryons, whereby we compute basic 
properties of holographic baryons, interaction with mesons, and ultimately nuclear 
forces. 
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What matrix? In our matrix model, the number of baryons is represented by the 
rank of the matrix. Therefore, for k-hodj baryons, it is a f/ (k) matrix model. If we 
path-integrate out the off-diagonal elements of matrices, we are left with k diagonal 
elements. It is these elements which represents the positions of the k baryons. In 
addition, there are a pair of complex k x Nf rectangular matrices whose classical 
values are related to the size of baryon. Together, they form the well-known Atiyah- 
Drinfeld-Hitchin-Manin (ADHM) matrix of instantons. 

Why matrix? In the large Nc QCD, mesons are open strings and light degrees of 
freedom, with mass ~ (9(1) while baryons are solitons with large mass ~ 0{Nc) 
. If we embed the large Nc gauge theory into string theory, baryons are described 
by string theory solitons, i.e., D-branes ^. Therefore, the dynamics of multi- 
baryons are described by a multi-D-brane system, which is nothing but the U{k) 
matrix model. Our matrix model is along the line of the ADHM construction [^, 
which is the matrix description for multi-instantons in gauge theory, or equivalently 
along the line of D0-D4 quantum mechanics. In the context of QCD, one must 
of course deform the description appropriately. Before we explain these changes in 
detail, let us first briefly review the holographic QCD. 

The application of gauge/gravity duality [|, ^ |10[ to large Nc QCD [|[], i.e., 
holographic QCD was developed very much in recent years, especially thanks to the 
D4-D8 model by Sakai and Sugimoto [O. This model starts with large Nc number 
of D4-branes compactified on a thermal circle []12|, representing pure QCD at the 
low energy, and incorporates Nf species of massless quarks by intersecting Nf pairs 



of D8- and anti-D8-branes |Tl|]. In the large Nc limit, the D4-branes are replaced by 
their dual geometry, and flavor D8 and anti-D8-branes are connected at the IR of 
the dual geometry, which is the geometrical realization of chiral symmetry breaking. 
The theory on the connected D8-anti-D8 brane in this dual geometry reproduces the 
low energy effective theory of light meson sector with remarkable accuracy. 

This theory has only two independent input parameters, and are free of am- 
biguities which low energy chiral Lagrangians such as Skyrme model |TB] possess. 
Furthermore, it contains not only the light vector meson but also infinite towers of 
massive vector meson, resulting in a model with infinite number of predictions. Most 
of these are irrelevant since at high energy the theory deviates from real QCD, but 
nevertheless there are many low energy processes which can be computed from this 
model. This aspect of the D4-D8 model [0, |1^ is one of the most important points 
that sets the D4-D8 model aside from the other more generic AdS-inspired models. 

Baryons in this D4-D8 setup turned out to be very interesting also. Baryons, in 
the large Nc limit, are inherently nonperturbative objects with their mass scaling as 
Nc [0 which is inverse of the genus expansion coupling constant. The holographic 
baryon in D4-D8 model is no exception, and can be introduced as instantonic soliton 

^The relevance of such a matrix model for baryons was previously emphasized by one of the 
authors i, Irf. 
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on the flavor D8-brane. Note that no new parameters are introduced in this step and 
all observables associated with baryons are computable by the D4-D8 model. ^ 

Following the original idea 0, ^, we recall that baryons can be also thought of 
as D4-branes wrapped on the compact S*^ that surrounds the QCD Nc D4-branes. 
Let us call them D4'-branes to distinguish them from those responsible for QCD. As 
such, its effective dynamics in the large limit involve open strings with both ends 
on the D4'-branes and also those connecting the D4'-branes and the D8-branes. The 
dynamics of k baryons would be U (k) gauge theory, and if we consider only the zero 
mode^ along the on which the D4' and the D8 are wrapped on, the theory on the 
D4'-branes reduces to a 0+1 dimensional matrix model. Motivated by this viewpoint, 
in this paper, we propose a new U{k) matrix model for holographic baryons. 

To illustrate our matrix model, let us remember a simple D-brane bound state 
of the Dp-D{p + 4) system and that we have two descriptions for this system. One 
is the solitonic description for the Dp-branes as instanton solitons in the D(p + 4)- 
brane gauge theory. This description is natural when gauge fields are weakly varying, 
namely at long distance scale with p ^ Is where p is the size of the instanton soliton 
and Is is the string length. The other is the open string matrix theory viewpoint, 
whereby the ADHM construction of the instanton soliton is naturally derived 



6[. This approach is more natural at a short distance scale, p ^ Is- For fully 
supersymmetric case, the two descriptions are equivalent for many purposes. 

Our matrix model for the k D4'-branes follows the latter viewpoint and is a 
D4'-D8 matrix theory compactified on a common S"^. There are several differences 
between our matrix model for baryons and the usual Dp-D(p + 4)-brane system. 
First, in our matrix model, there is a Chern-Simons (CS) term (supersymmetric 
versions were studied in [pTj This term originates from the fact that there is a 
Ramond-Ramond (RR) flux on the on which the D4'-branes are wrapped, and is 
in fact the same type of term that allowed Witten to identify wrapped D5-branes as 
baryon vertices in AdS/CFT description of maximally supersymmetric Yang-Mills 
theory 0. For us, this term turns out to play a crucial role in dynamics of baryons 
in general and in the baryon-baryon interactions in particular. Second, since the 
D4'-branes are living on the D8-brane at the IR bottom of the warped geometry, we 
need to take into account the warped geometry to derive the U{k) matrix model. The 
point is that if we put the k baryons at short distances, the warp factor approaches 
almost constant values, therefore the effects of the warped geometry is simply just 

^Employing this picture, static properties |l^, interactions with mesons [15[ |l^, electromag- 
netic form factors and more recently nucleon-nucleon potential |2^, ^| have been derived 
again with remarkable accuracy. Although most literatures considered nucleons, higher isospin 
baryons can also be treated on equal footing |2^, . 

^The subtlety of taking of a decoupling limit of the non-zero modes along is a long-standing 
problem in holographic QCD. As a result, it is obscure why in low energy effective theories such 
heavy objects like baryons are well-described. For one possible explanation concerning supersym- 
metries offered by one of the authors, see p3. 
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rescaling of a coupling constant plus mass terms. This in turns implies that the 
supersymmetry is broken explicitly by the dual geometry |T^, ^ of the QCD D4- 
branes. In practice, the supersymmetry breaking would manifest as various mass 
terms and potential energy in the would-be vacuum moduli space. 

One motivation for us is that such a matrix formulation can easily accommodate 
a large number of baryons. One of more grey area of nuclear physics, from the purely 
theoretical viewpoint, is how one handles many-nucleon systems such as nuclei. The 
fundamental theory of QCD is even less effective there. The matrix model of ours 
can be written down immediately for all k, and should contain in principle, not only 
2-body interactions but arbitrary fc-body interactions built-in from the beginning. 
This may open up a novel window for dealing with finite k nuclei and large k physics 
such as necessary for neutron star, in term of strongly-coupled QCD. 

Another motivation, which is perhaps more practical for now, is already apparent 
in the holographic soliton picture of baryon |jl5|, |l^ . Remember that the holographic 
size of the solitonic baryon is ~ 9.6/(MkkV^), where A is the 't Hooft coupling and 
Mkk ~ 1 GeV is the mass scale of the vector mesons |jT5|, |16|. Although this size 
is larger than the Compton size of the baryon ~ 1/(MkkA'^cA), obeying the usual 
criterion for validity of field theory soliton picture, it is unfortunately as small as 
the local string scale if^ ~ 1/(Mkk"\/A)- One should be dubious whether the field 
theory soliton is really accurate enough for such a small object. The supersymmetry 
breaking scale is ~ 1/Mkk, so the supersymmetry is approximately valid at the 
very short distance such as the string scale, so one hope for relative stability of the 
baryon physics in interpolating between the string scale and the QCD meson scale 
p^ . Nevertheless, we cannot expect quantitative agreement for all observables. 

Open string description of D-brane interactions is deemed to be relevant for dis- 
tance scale below string scale, so we have a motivation and an opportunity here to 
reexamine holographic baryons from a D4'-D8 (compactified on 5*^) matrix descrip- 
tion. With supersymmetry broken at M^k and below, one should in general expect 
the solitonic viewpoint and the matrix viewpoint would disagree quantitatively; it 
would be interesting to compare our matrix description with the previous soliton 
description. 

In section 2, we drive the matrix model by starting with D4'-D8 gauge theory 
compactified on S"^, and extracting leading order supersymmetry-breaking effects at 
very short distance. The action we find is valid at distance If^ and below. To show 
the effectiveness of our matrix model, we display only two simple examples in this 
paper: energy functions of static configurations for k = 1 and k = 2, for one and two 
flavor (s). 

In section 3, we discuss k = 1 and determine the holographic "size" of the 
baryon, value of which affects numerous observable quantities. This parallels the 
energy function estimates for the solitonic baryon in spirit but differs in detail. For 
two or more flavors, we find the energy function with a different numerical coefficient. 
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leading to a new estimate for the holographic size which is larger by a factor of 
(5/4)^/"^ than the soliton result. We discuss its implications. With a singe flavor, 
for which the soliton model has no computation due to lack of finite and smooth 
self-dual U{1) instanton solutions, we find a smaller but still non-vanishing size. 

In section 4, where k = 2 is studied, we compute a baryon-baryon potential 
at short distance for two flavors. Interestingly, integrating out the auxiliary gauge 
potential in 0+1 dimension turns out to give a universal repulsive core of the nuclear 
force. The core consists of three terms; the isospin-independent central term, isospin- 
dependent central term, and isospin-dependent tensor term. The computation here 
is far simpler than the soliton computation, yet gives the same type of results, except 
that numerical coefficient of the latter two turn out to be larger by a factor of 5/4. 
We explain this in simple terms based on the above k = 1 results. We also show that 
the repulsive core is universal for any baryon state for the two-flavor case. 

In section 5, we close with discussions on possible ramiflcations of this new model 
of baryons and baryon dynamics. 



2. The matrix model 



We derive the action of the matrix quantum mechanics in this section, which will be 
used for getting the baryon spectrum (section the nucleon-nucleon potential and 
the universal repulsive core (section ^ . 

The matrix model is nothing but the low energy effective action on the k D4'- 
branes embedded in flavor Nf D8-branes in Witten's geometry. Those who are not 
interested in the string-theoretical derivation are advised to see only the action (in 
section |2Tl| ), skip the rest and go directly to section ^ 

2.1 Action 

The matrix model action we derive is a U{k) quantum mechanics. 



S 
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dt tik 



(DoX 
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4A2M^K 
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D-T'^^X^'^Xaa + D-T 



+N, / dt trfcAo , 



(2.1) 



Here A = NcgqQD is the 'tHooft coupling constant, and Mkk is the unique dimension- 
ful constant. The dynamical flelds are X^ and w, while D and Aq are auxiliary flelds. 
All the flelds are bosonic. We claim that this matrix model describes the fc-baryon 
system, according to the holographic principle in string theory. 

Our 1-dimensional matrix model is a deformed ADHM matrix model. The 
ADHM matrix model has been extensively studied in the context of D-branes in 
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Figure 1: Fields in the matrix model. 



string theory and instanton calculus (for a concise review, see |^). Our theory is 
deformed in the following two points: 



• Addition of the CS coupling. The last term of the action (|2.1|) is a CS term in 
1 dimension. 

• Mass deformation. A part of the dynamical fields, and w, are massive, in 
contrast to the standard ADHM matrix model. 

We will describe how the ADHM matrix model and these deformations appear in 
the holographic QCD, together with the detailed derivation of the coefficients in the 
action. Note that in the absence of the field w, the CS term and the mass term, our 



matrix model looks close to the BFSS Matrix theory for M-theory ||3l[ or the IKKT 
matrix model as the integration of the auxiliary field D results in a potential of 
a commutator type, tr([X, X]^). 

The symmetry of this matrix quantum mechanics is 

U{k) X SU{Nf) X S0{3) (2.2) 

where the first U (k) is a local symmetry with which the gauge field Aq is associated, 
and the remaining SU{Nf) x 5*0(3) is a global symmetry, k refers to the number 
of baryons of the system, Nf is the number of fiavors in QCD, and SO {3) is the 
rotational symmetry of our space in which the baryons live. The notation of the 
action is better-understood if we embed the rotational symmetry as 

^0(3) C S0(4) ~ SU{2) X SU{2) (2.3) 

where the additional dimension in fact corresponds to the holographic dimension. 
This S0{4) is broken down to the S0{3) by the mass deformation. 

We summarize the representation of the fields in the table 1. The indices with 
respect to the U{k) gauge group are implicit. In the action, the trace is over these 
U{k) indices. One can think of the fiavor symmetry as U{Nf), while the overall U{1) 
part of it is identical to the overall U{1) part of the gauge symmetry U{k), as seen 
in how they act on the bi- fundamental field w. 
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In the action, the definition of the covariant derivatives is DqX^ = OqX^ — 
i[Ao,X*^], Dqw = OqW — iwAq, Dqw = OqW + iAqw, and (s = 1,2,3) is the 
Pauh matrix. The spinor indices of X are defined as Xaa = X^\aM)aa and = 



X (aM)""^ where = (^^, 1) and au = {—i^, !)• We follow the notation of [BO 



2.2 Derivation in gauge/gravity duality 



Our matrix model (|2.1|) is nothing but a low energy effective field theory on D-branes. 
The D-branes of our concern are D4'-branes wrapping 5*^ of a background geometry 
given by Witten |T2| (the metric originally given in [^). This D4'-brane is called 
"baryon vertex" |^ Q in the gravity side of the AdS/CFT duality, which corresponds, 
as the name shows, to a baryon in the field theory side. As we are dealing with a 
fc-baryon system, we place these k D4'-branes close to each other. We are going to 
derive the effective action of this collection of the D4'-branes wrapping the 5*^, via a 
standard technique in string theory. 

The D-brane action is affected not only by the geometry and the background flux, 
but also by the presence of the probe Nf D8-branes which are responsible for quarks 
a la Sakai and Sugimoto . We have additional strings connecting the baryon D4'- 



brane and the flavor D8-brane. From the viewpoint of the baryon D4'-brane effective 
field theory, this string provides the field w in the bi-fundamental representation. 
Together with the field X'^^ in the adjoint representation whose diagonal eigenvalues 
specify the location of the D4'-branes in the transverse directions (but longitudinal 
to the D8-brane worldvolume) and thus the location of the baryons in our real space, 
the matrix model action is written. 

Since the D4'-branes are completely inside the worldvolume of the fiavor D8- 
branes, the matrix model is very close to the so-called supersymmetric ADHM matrix 
model which is nothing but the effective action of k D-instantons on Nf D3-branes in 
flat spacetime. The deformation is due to the curved geometry and the flux, which 
break the supersymmetry explicitly. Here we keep only the bosonic flelds, which is 
enough for computing classical quantities. 

Although the D4'-branes (and the D8-branes) wrap the 5*^, we trivially reduce 
the dependence (a dimensional reduction with the assumption of no dependence 
along S^), so that the resulting action is in one dimension, i.e. only time direction. 

The background geometry and the flux given by Witten [0 (and Gibbons and 
Maeda [^) are written as 

ds^ = {U/Rf/\r]^,dx^dx'' + f{U)dT^) + {R/Uf^^{f{U)-^dU^ + U^d^l) ,(2.4) 
e^ = g,{U/Rf\ F, = dCs = ^e,, (2.5) 

where f{U) = 1 - U^JU^, and = ngsNJl I/4 = 87rV3 is the volume of the S\ 
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and €4 is the volume form on it.*^ The r direction is compactified with the period 
r ~ r + 27r/MKK where Mkk = (3/2)?7j^'^i?~^''^, so that the geometry is everywhere 
smooth. The relations to the QCD variables are 

R' = l'-^ll U^K = UuNMl 9s = ^^^. (2.6) 
Z Mkk y 27r Mkk h 



Convenient coordinates used in O, in 



are 



U' = U^^ + UKKr', O^l^r, y + ^ re^^ (2.7) 

The flavor D8-branes are located at y = 0. 

Below we derive the action, by looking at, flrst, the background RR flux, and 
second, the effect of the background geometry. 

Chern-Simons term 

The important term is the last term of the matrix model (|2.1| ), which is a CS 
term in 0+1 dimension. Using the background RR flux ( p. 51) ^ we can compute the 
CS term on the D4'-brane as 

^ = ^2^ / "^'^ ^'''"''"^C,,,,,,F^p = i/"^^ / 

= nJ dt tiAo . (2.8) 



Note that the overall factor A^^^ shows that the D4'-brane should be supplied with A^^ 
fundamental strings (the end point of the D4'-D8 string serves as an electric charge 
on the D4'-brane), which means that the D4'-brane is indeed a baryon P, 

Interestingly, this 0+1-dimensional CS term was used in ||2^ and for the 



ADHM matrix model for supersymmetric dyonic instantons with a CS term in 5 
dimensions. There, the CS term was argued for heuristically (p^ studied fermions 



and anomalies to reach the 0+1-dimensional CS term). Here we have derived the 
0+1-dimensional CS term from string theory. 

Mass terms and overall normalization 

The Dirac-Born-Infeld part of the action for a single D4'-brane is 

S = -Td4 j <fi e-^ V-det(GMJV + 27ra'FM7v). (2.9) 



^We use the standard normalization for the forms, C3 = (1/3!) Ctjkdx' A dx^ A dx^ and — 
(1/3!) dcC^jkdx°' A dx' A dx^ A dx'' . 

^We are using the normalization of the RR field in which the RR charge is measured in units of 
27r, see Appendix A of [|ll[ . 

^This Chern-Simons term is also important in producing the correct statistics of baryons, which 
should be either fermionic or bosonic depending on whether is odd or even. For detail of nucleon 
statistics in our matrix model, see Ref. IBSl. 
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We consider a D4'-brane situated at ?/ = which wraps the S*^. Then 

S=-^ fdtV, {{R/Uf'U'Y {U/nr'/'y^o (2.10) 
9 s J 

with V4 = Stt^/S the volume of a unit four-sphere, where the induced metric is 

= - {tT " ^ 5 iff' ■ P.ii) 

The index of X* runs for our 3-dimensional space, i = 1,2, 3. So, we obtain 



s = -^jdt ^^R'u^i - idox^y - l^idozy . (2.12) 

with Tdp = 27r/(27rv^)P+i. 

We expand this for small Z and small X. Using the expansion U = Uk_k{1 + 
{l/3)U^lZ^ + 0{Z^)), and a redefinition 



2 f R 



3/2 

X^ = ZI^^\ z, (2.13) 
we obtain a quadratic Lagrangian 

1 . . 1 . . 1 . .J 

(2.14) 



21 TX 



-1 + \{d,X^f + ^(5oX^)' - ^M^k(X^)^ 



In this way, the overall normalization of the matrix model action,'' as well as the 
mass term for the field X^, are provided in the matrix model ( ^.ip . This mass term 
originates entirely from the expansion of U sitting in front of the square-root in 

(in- 

This leaves the mass term for w. We assume that it is 

AA^eMKK /■ , 1 



547r 



jdt^Ml^w^w^, (2.16) 



to be added to ( p.l4|) , which completes the quadratic part of the matrix action in 
( p.l|) . This (|2.16|) is an educated guess via a comparison with the sohton picture 



^Note that higher order terms in X"^ are not suppressed by 1/A. In fact, the leading correction 
to the mass term is 

+ \Ml^{X^f-\MiAX^f + --- (2.15) 

In the foUowing, we assume that the magnitude of X"^ is smah so that we can ignore the higher 
corrections. In fact, as we consider the wave function of this X'^ after the quantization around the 
vacuum of the matrix model, we obtain a Gaussian wave function with the width suppressed by 
1/Nc which is quite small, and this approximation is valid. 
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(though in principle one can derive this by computing string scattering amphtudes). 
In fact, as we will see in the next section, the mass term in ( p.l4| ) coincides with 



the soliton picture, and the following argument suggests that w mass also coincides. 

Recall that (X*,X^,W7) form the ADHM data for instantons. When the instan- 
ton size is very small compared to the supersymmetry-breaking scale 1/Mkk, the 
background geometry is effectively flat with approximate supersymmetry. The D4'- 
D8 system will inherit supersymmetry, and thus the data (X*, X"^, w) is equivalent to 
the supersymmetric instanton on the D8-branes; At this zero-th order, the instanton 
picture of the D4'-brane on the D8-branes is still valid. Note that this is not yet 
the baryon but merely a purely magnetic instanton solution. Then, we turn to the 
effect of supersymmetry breaking at Mkk scale, and evaluate the potential energy 
of the instanton at the first order. This is in fact one way of obtaining the mass term 
for X*^ in (|2.1| ). Similarly, using the fact that \w\'^/2 is the size squared for a single 
instanton, we find ( p.l6| ) for the matrix model.^ 

Commutator term 

Finally, let us compute the coefficient in front of which is related to the 
coefficient in front of the famous commutator term [X, X]^ in D-brane quantum 
mechanics. 

Let us expand the generic Dp-brane action {\2.9\) to the quadratic order. 



S = -T^pj<F^^i e-'t'^/^d^iG^\{27:a'fFMNFpQG^'''G'''^ . (2.17) 

We now make a dimensional reduction to get the commutator term from the YM 
kinetic action. The relevant formula for the dimensional reduction is 2TTa'AM = 
X^Gmn for diagonal metrics. Then, the action is 

L oc 2G^^G^jDoX'DoX^ + 2G°^G,,DoZDoZ 

-[X\ X^] [X^ X']G..G,,^^ - 2[X^ Z] [X^ 



(X (DoX0^ + (Z^oX^)2 + i^^([X\X^]2 + 2[X\Xt) 



R 



= iDoX'r + -^^X'M^^[X^^,X'^]\ (2.18) 
On the other hand, the commutator term can be written by using the auxiliary 



*One may recall that the soliton representation of the holographic baryon has a Coulombic 
potential as well [ p^ [l^ . This is not apparent yet in the matrix model but will arise by integrating 
out the gauge potential Aq of D4'. The crucial difference from the mass terms is that the Coulombic 
energy arises from a quadrature of an electric field (in the soliton picture) which is itself a first 
order deviation. We will see later that, perhaps because of this, the Coulombic energy differs 
quantitatively in the two pictures. 
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field D. If we start from the action of 



S = c dt tT 



(2.19) 



then by integrating out the field D we obtain 

1 



S = c dt ti 



167r2a'2 



(2.20) 



Comparing this with the normalization we obtained in ( p^.l8[ ), we obtain the expres- 
sion for the commutator term of our matrix model as 



XN^Mkk 

54:71 



dt tik 



(2.21) 



The w term coupled to D is written down in and we can just use it with the 
same normalization as the DXX term. 



3. Single baryon 



In this section, we study the k = 1 case, i. e. a single baryon. Our quantum mechanics 
directly gives a spectrum of the baryon. 

First, we evaluate the Hamiltonian of the quantum mechanics with k = 1. Then 
we analyze the vacuum of the system for Nf = 1 and Nf = 2 respectively. For Nf = 2 
system, we calculate the baryon spectrum. The computation of the spectrum, as well 



as the quantization procedure, closely follow the soliton approach of |jT5|, |T6|, although 
the derivation of the Hamiltonian is different. Finally we discuss meson couplings. 



3.1 Hamiltonian 

Let us compute the Hamiltonian for a single baryon k = 1, with generic Nf. 

First we explicitly integrate out auxiliary fields D and Aq. As for the terms 
including the field D, since for k = 1 the field X is now not a matrix but a number, 
all the X couplings drop off, and we obtain 



S 



XNrM, 



KK 



D 



547r 
XN.Mkk 



547r 



36772 
36772 



+ (./;i(K)*)' + (^(^)*)'-2K(K)XK)1 , (3.1) 



where, the first sum is over three Pauli matrices and we omit the dots in the dot- 
ted spinor, a = 1,2. 5*^ gives a so-called ADHM potential. Minimization of the 



^In teol the field D is anti-Hermitian, while our D is defined to be Hermitian. 
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ADHM potential is equivalent to the ADHM constraint, which should be solved for 
construction of instantons in the ADHM formalism. Note that since Lagrangian is 
Hermitian, f"'^X^°'Xaa and T^^w^Wai are real. 

Since our theory is 0+1-dimensional, the gauge field Ao is an auxiliary field, and 
we integrate it out explicitly. The terms including in the matrix model action is 



XN.Mkk , 



5Att 



dow'^{-i)waAo + iAoiv'^doWa + {AoYiv'^Wa + ttt— 

AMkk 



(3.2) 



So, the equation of motion for this Aq, in other words, the Gauss law constraint, is 



547r 



AM] 



KK 



(3.3) 



Then after path-integration over ^40, we obtain 



Awfwl V AMkk 



547r 



■ (3.4) 



Using the definition of the momentum conjugate to the field w 



dWn 



547r 



<9o< 



Awjwi^ \AMkk 



547r 



+i{wid,wl-d,wiwl)]iwt 



(3.5) 



we obtain the Hamiltonian 



H = 



547r 

+^;^ [HiwiTwiiwir + {w{{w{rf + {wHwirY - 2^(k)x(^^)*] 



3%2 



547r 



+ (<^o<-^o<0' 



(3.6) 



3.2 Single flavor 



Let us minimize the Hamiltonian to find a vacuum of the k — 1 system. We consider 
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first the case of tfie single flavor, Nf = 1. We put the following ansatz,^° 

Wa=l = Pi , Wa=2 = P2 , (3 

where p's are real constants. Then, the Hamiltonian is 

XN.Mkk 



H 



547r 



1 f 27n \' 2 1 ..2 2 4A2M^K 4 



2 V AMkk/ ^ 3 ' 3677 



(3.9) 



where we have defined 2p^ = p\ + p^. 

Each term in the Hamiltonian (|3.9|) has a physical meaning. 

• The first term oc p~^ is induced by the CS term and the Aq path-integration, 
and it can be interpreted as a self-repulsion of the dyonic instanton in the 
soliton picture. As described in |]15|, |16|, the instanton has an electric charge, 
so the self-energy should be lowered by expanding the size of the instanton, 
thus resulting in a negative power in p. 

• The second term is from the mass term of our matrix model, thus comes from 
the curved spacetime of the background. In terms of the instanton, the location 
of the instanton along the direction z affects the total mass of the instanton, 
due to the curved geometry. 

• The third term oc p^ is from a path-integration over auxiliary D fields which 
for example, yielded the commutator square term in the matrix model action. 
So this corresponds to the ADHM potential term. For the single flavor, there 
is no U(l) instanton except for the small instanton singularity in fiat space, 
and this term ensures it, in the absence of the dyonic coupling and the curved 
geometry. 

The Hamiltonian is minimized at a nonzero p, but the minimization problem is 
a non-linear equation. With a help of the fact that we are working in the large A 
limit, we reduce the problem to a linear one. We put 

p = xA°Mkk (3.10) 



^^A more general ansatz is, 

= pie^(-i*+^i) , «;^=2 = P2e^'^^*+^=) , (3.7) 

where p's, w's and s's are real constants. However, this results in the vacuum which is the same as 
what we will find below in this section. In fact, if we minimize the Hamiltonian H with respect to 
vi and V2, we find vi = V2, and the Hamiltonian H is independent of v. This is interpreted as a 
manifestation of the gauge invariance. As a result, we can always choose a gauge vi — V2 — 0. Note 



that if one choose another gauge, for example, ^0 = 0, then the Gauss law constraint (3.3) forces 
w to have a time-dependence, i.e. wi = W2 7^ to satisfy the Gauss law constraint (3.3). Simplified 



ansatz (3.8) is not consistent with a generic gauge choice for Aq. This is in good contrast to the 
situation of the ADHM vacuum for supersymmetric Yang-Mills-Chern-Simons instanton studied in 

MM- 
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where a; is a constant coefficient, and a is a constant power. Then, each term in the 
Hamiltonian ( |3.9| ) scales for large A as 

1 / 277r \ 2 x-2-2a ^ A/T'^ ^2 x2o ^A^M^ . ,4a+2 /"i i i \ 
p ~ A , -MkkP ~ A , p ~ A . (3.11) 



2 VAMkk; ' 3 ' 3671 

In minimizing the Hamiltonian, the ffist term, which has a negative power of p, should 
be balanced with either the second or the third term which has a positive power of p. 
If it is with the second term, then from the above A-scaling we need to have —2 — 2a = 
2a, thus a = —1/2. However this value means that the third term has a power larger 
than the ffist and the second terms, so leading to an inconsistency. Therefore, we 
conclude that the minimized Hamiltonian is dominated by a cancellation of the ffist 
and the third term. This leads to —2 — 2q; = 4q; + 2 which is solved as a = —2/3 at 
which indeed the second term gives a smaller contribution, which is consistent. For 
this reason, we can safely ignore the second term at the large A, so 



XN.Mkk 

H ^ 

547r 



1 / 27n y 4A^M^K ^ 

2 VAMkk; 3H 



(3.12) 



Ignoring the second term means that the ADHM-like potential (the third term) is 
much larger than the curvature scale of the background spacetime. This is natural, 
since the ADHM-like potential has the scale of the string length, in D-brane effective 
actions. We will see in the next that for the two flavor case the ADHM-like potential 
can vanish so that finally the Hamiltonian is minimized by the cancellation of the 
ffist and the second terms. 

The value of x minimizing this Hamiltonian is computed, as 

p = 2-2/39y^A-'/=^MK^ . (3. 13) 

The minimized value of the Hamiltonian is 

H^i^ = 2''/^X'/'N,MKK- (3.14) 

3.3 Two flavors 

Let us consider the more realistic two-flavor case. First, to eliminate the contribution 
from the ADHM potential term (the third term), we need to satisfy the ADHM 
constraints, r'^^w^Wai = for all Pauli matrix directions, or equivalently, 

Nf Nf Nf Nf 

E<=i«=2)* = E<=2«=i)* = 0' Ei<=ii' = Ei<=2i'- (3-15) 

i=l i=l i=l i=l 

Once this condition is met, the ADHM potential disappears, and the total energy is 
lowered drastically as the power in A changes. This can be achieved by the following 
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generic choice 



^r = {^ , ^T' = (^'] . (3.16) 



0/ . \P 

'a ^ ' ^ 

Note that this is a generic solution minimizing the ADHM potential, since the con- 
dition is invariant under the U (2) global transformation on the spinor index a and 
the t/(2) flavor symmetry, 

<->^.v«[f/;]/- (3.17) 



Then the Hamiltonian, after we include X^-dependence as well, is 
AAT.Mkk 



E 



(3.18) 



547r 

This is minimized at 

p = 2-1/^3^/^ . (3.19) 

The minimized value of the Hamiltonian is 

H^in = Q-'^^N,MKK. (3.20) 

This is independent of A, thus in the large A limit, we see that -ffmin for the two flavor 
case is far smaller than that of the single-flavor case. Recall that the classical mass 
of the holographic baryon is AA'cMKK/277r + H^am where the first term comes from 
the constant part of ( p.l4| ). 

Here, the variable p is nothing but the instanton size in the soliton approach, 
since we have chosen a correct normalization for this p, a la ADHM formalism in the 
flat spacetime. Let us compare our Hamiltonian with the one obtained in the soliton 
approach, |15|, |l6l, where a potential for moduli of the single instanton solution 
was computed. The Hamiltonian |l6l for the instanton size modulus p and the 
instanton location Z along the direction is 



-'^soliton _ . 



9 ■ 3^7r^ 11 9 



L5(AMkk)'p' ■ 3""""'" ■ 3" 



(3.21) 



again without the rest mass term \NcMky^/2'Jtt. We first note that the quadratic 
terms in and in p coincide with ours. The coincidence of the mass term is 
nontrivial, while the one for p is not accidental, since we have computed the p mass 
term in ( |2.1| ) by resorting to the soliton picture in the supersymmetric limit. In some 
sense, they are first order terms whose evaluation used the zero-th order solution. 
The real comparison is with the first terms which are proportional to 1/A^ (times the 
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naive baryon mass computed from DA' tension); the two are structurally identical 
but the matrix model result is larger than the soliton result by a factor of 5/4. 

In the soliton picture, this term arises as the five- dimensional Coulomb energy 
associated with U (1) baryon charge; the latter is, due to the holographic map, a gauge 
charge of real gauge field on D8. In the matrix model, it comes from integrating out 
the (non-dynamical) Aq gauge field on DA'. Either way, it comes from a quadrature 
of an excitation field of order 1/A, which suggests that, in the supersymmetric limit 
of very small instanton, the Coulombic energy captures the deviation from the zero- 
th order "instanton = ADHM" configuration more effectively. We suspect that this 
explains the numerical difference. At any rate, 5/4 is fairly close to 1, implying that 
the difference between the two approaches is relatively minor. Since the two such 
models in general live in two vastly different validity regions, respectively, one could 
have expected a larger difference. It is the presence of approximate supersymmetry in 
distance scale from 1/Mkk down to If^ and below that give us this relative stability. 
The baryon physics of one does not deviate a lot from the other. 



3.4 Quantization 

With this result on the matrix model vacuum in mind, we can now quantize the small 
fluctuations of the k = 1 Nf = 2 matrix model. This spectrum should correspond to 
the baryon spectrum. 



As we have seen, what is different from |]T6| is just the coefficient in front of the 
1/p^ term in the Hamiltonian. So, we can just track the difference in the computa- 
tions of and find the spectrum of our matrix model. This difference reflects in 
the constant Q in ||16|, which is now multiplied by 5/4 in our case. Then the mass 



formula for the baryon excitation is 



M = Mo + ^ + + ^ ' . (3.22) 

Here Hp, = 0, 1, 2, ■ ■ ■ and I = 1/2 = J/2 with spin J and isospin I. The difference 
from |T6| is just the coefficient of A*"^ in this expression. Mq is the mass of the D4'- 



brane which is equal to the first term in ( [^.14D , AA^c^kk/STvt. 



This formula in particular means that we can obtain 

M;=3 - M;=i = 0.5693Mkk • (3.23) 

If we use Mkk = 945 [Me V] which is fit by the p meson mass, we obtain 

Ma - Mn/p = 540 [Me V] . (3.24) 

This is larger than the experimental value 292 [MeV]. The situation is similar to the 
soliton approach [|l^ which gave the value 569 [MeV]. 
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3.5 Meson couplings 

As we noted toward the end of subsection (|373|), tlie energetics of p liere differ 
quantitatively from the previous estimate based on the instanton soliton picture 
in D8-brane gauge theory. One consequence is that the classical value of p, which is 
the holographic size of the nucleon, is slightly larger that its previous estimates in 
Refs. JlS], 0. Denoting the latter by Psouton, we found 

= y|psoiiton- (3.25) 

This new estimate modified the baryon spectra as we just saw above, but it should 
also affect couplings to mesons. 



This can be seen most clearly from the effective action approach |T5|, |17|, where 
a tree-level effective action in the five- dimensional bulk captures nucleon M coupled 
to mesons 



J 3e^ 



2^2 

(3.26) 



rrij^ and are known functions of the holographic Z-coordinate, determined by the 
dual geometry. Actual four- dimensional nucleon is the lowest lying mode of five- 
dimensional Dirac field (here, denoted by the same symbol in abuse of notation), 
while the infinite tower of mesons are embedded into the five-dimensional flavor 
gauge field ^^^^^ = A^'^^^ + A^^'^'^^ as in Ref. [O. Note that the size parameter 



Psoiiton appears explicitly only once, in the last term, so the holographic size of the 
nucleon will affect couplings to SU (2) iso-triplet meson coupling with specific chiral 
or tensor structures only. 

Previous estimates of the meson-nucleon-nucleon couplings were based on the 
instanton soliton viewpoint and may not be completely compatible with our new 
matrix model. Nevertheless, we note that these couplings were read-off entirely from 
the long distance gauge field configuration associated with the baryon |jl5|, |^ , which 
has to he the case since the Compton wavelength of mesons at ~ 1/Mkk o'^e much 
larger than the size of the soliton core Psoiiton ~ 1/MkkVX- This suggests that the 
derivations of meson-baryon coupling are generally safe from short-distance physics. 
The one place where this reasoning can go wrong is the size of the soliton itself, 
whose estimate relies heavily on short distance physics and which in turn affect the 
long distance magnetic field of the soliton. Thus, we expect that the effect of the 
latter manifests in the meson-baryon couplings via the single quantity, p^. 

We can classify the cubic meson-nucleon-nucleon couplings from (p.26|) into two 



classes. The first class consist of those whose leading large Nc behaviors originate 
from the minimal coupling to A'^'''^\ This includes dimensions-four couplings to 
vector mesons, and dimension-four couplings to all iso-singlet mesons: vector, axial 
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vector, or pseudo-scalar. With the exception of those to iso-singlet vectors, this class 
of couplings are subleading in l/N^. Let us collectively denote them as g'^^\ 

The second class, g'^^^\ have the leading Nc contributions arising from the last 
term of ( p.26| ) and thus proportional to p^. This class consists of coupling to iso-triplet 
Goldstone boson (namely pions) , minimal couplings to all iso-triplet axial vector 



mesons 



p!7[] , and dimension- five tensor couplings to all iso-triplet vector mesons . 



When we replace Psouton by p^, we therefore find 



9 — fl'soliton' 9 — \/ ^ 5'soliton • l-J-^'J 

Numerically, the latter represents about 12% increase for g^^^^^s. 

Translating to quantities more directly related to data, we have for example the 
charge form factor and the isospin-independent part of nucleon-nucleon repulsive 
core unaffected, whereas the magnetic form factor (thus the anomalous magnetic 
moments also) and the isospin-dependent part of nucleon-nucleon potential would 
be increased by factor of 5/4 in the large Nc limit. See the next section for a direct 
matrix model computation of the nucleon-nucleon potential at short distance, which 
confirms this expectation. 

Whether or not such a shift is beneficial in reproducing QCD with iV^ = 3 is 
unclear and needs to be studied more. We would know the answer only after we have 
classified and computed subleading corrections in this D4-D8 model. Regardless, we 
note that, purely within the context of studying baryons in the D4-D8 model, this 
represents a significant change to the leading l/Nc computation (and not a subleading 
correction) which could have not been obtained by any other method we know of. 



4. Two-body baryon interaction 

The baryon interaction at short distance can be obtained by classically integrating 
out Aq, as in the case of the single baryon. Now, with two baryons, we have the 
matrices charged under the U{2), so, in particular, Aq has four components, 

Ao = A°i2x2 + Ay + Ay + Ay . (4.i) 

Recall that, in the previous section, solving for the equation of motion of the overall 
f/(l) gauge field gave the 1/p^ potential for individual baryon. Now we are interested 
in integrating non-Abelian part of Aq as well, which should generate interaction 
energy for a pair of baryons. 

Note that Aq {a = 1,2, 3) does not show up in our CS term in ( |2.1| ), because the 
latter contains tiAo only. Terms including Aq appear in the kinetic term of X and 
that of w. 
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4.1 Two-baryon configuration 

In order to evaluate the action, first we fix tlie vacuum of tlie ADHM potential in 

(O- 

4.1.1 Single fiavor 

The classical single baryon configuration is specified by the vacuum configuration of 
w. In section |3.2| , we obtained the classical vacuum of the matrix model as 

Wo. = u(\\p, (4.2) 



/. 

where ?7 is a 2 x 2 unitary matrix. This representation is equivalent to having 
Pi, p2 in the notation of section |3.2| . This unitary matrix comes from the global 
symmetry ( p. 17] ). Now we have two baryons, so the baryons are specified by two 
unitary matrices, 

<^^ = f/^^^(j)p, wr = U^''(\)p. (4.3) 

One can work out the baryon interaction potential with this, but we will perform 
the computation only in the 2-flavor case in the following, since it is realistic. 

4.1.2 Two fiavors 

For the two-flavor case, the vacuum configuration of the potential is given by the 
integration of the [/(A;)-adjoint field Dab in the matrix model (|2.1|) , which is the 
ADHM constraints, 

f (X'^^X,^ + w^^w^^ = . (4.4) 

Here we explicitly write the baryon index A,B = l,...,k. In this section we treat 
two instantons so k = 2. The generic ADHM configuration satisfying this equation 
is nothing but the ADHM data of two YM instantons. It is given by 

w^r' = Uif='^ Pu w^r' = Uif='^ p2 , (4.6) 

where the locations of the two baryons are given by the diagonal entries in Xm so 
that tm is the inter-baryon distance, and 

Ym ^ -^tr [cTMrNcrM {{U'^'W^'^ - (f/(^))tf/(i))] . (4.7) 

Here U^^^ and f/^^) are SU{2) matrices which denote the moduli parameters of each 
baryon, and ctm = {ir, 1), ctjv/ = {—ir, 1). 
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This ADHM data was explicitly used in the soliton approach |^0|]. In terms of 
the YM instanton, these degrees of freedom are gauge rotations of the "flavor" gauge 
group, and after the quantization, they become the spin and the isospin of each 
baryon. They can be written by real unit vectors a^^ and 



(1)- 



2x2 - 



(2)- 



2x2 ■ 



(4. 



with (ai'V + (a?V 



.(2)a2 



1. This is the correspondence to the 
notation of [|TB]. Using this, we obtain expressions which will be useful later, 

rMYM = 0, (4.10) 



pIpI 



■tr 



(f/(i)tf/{2)_f/(2)t^(i)y 



pIpI 



1 - io^o^f) , 



tr [t/(i)tf/(2)] = tr [t/(2)tf/(i)] 
4.2 Baryon interaction potential 



(4.11) 
(4.12) 



We shall integrate out the matrix and compute the two-baryon interaction po- 
tential. For integrating out the SU{2) components of Aq, we first write down all the 
terms including those components in the matrix model action ( p.l|) . First, 

triDoX'r = 2 {{Airrl + (A^)2(r^, + AY^,) + 4(^^)2^,^) - SAIAItmYm (4.13) 



Here we used only ( [4.5|) . But if we further use the explicit expression (^T^) and the 
vacuum for w ( [4.6|) , we obtain ([4.10|) , so the last term in (|4.13| ) vanishes. Next, the 
w kinetic term is 



trDo</^o^„. = 2(p? + pI) ((^0)2 + (Air + (Air + {Alf) 
+4p,p,KAl tr [t/(^)t^(^)] + 4(p? - pl)A',Al 



(4.14) 



In these kinetic terms, the component Aq appears only as a form (Aq)^- Thus we 
can minimize it independently with A^ = 0, meaning that we can just ignore the 
component Aq. So, the total kinetic action plus the Chern-Simons term is 

XN.Mkk 



5Att 

AAT.Mkk 
54-11 



j dt tr [{DoX^^f + DowfDow^i] + j dt tiAo 

2(AJ)M. + SiAlfY!, + 2{pl + pI) {{Alf + (Air + (AlY) 



dt 



lOSvr 



+4p,p,A',Al tr [f/(i)tt/(2)] + A{pl - pl)A',Al + ^A\ 



KK 



(4.15) 



""^""^To obtain the first equality of (4.11), we have used the following formula for unifying the double 
trace: 

(tr [ijMibo + ihT')]y = 2 tr [(6° + tbW') (6° ~ tb^r^)] . (4.9) 
Then, using + (af ^)2 = 1 and (of ^)2 + (af = 1, the second equality of follows. 
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As the action is quadratic in the remaining components of Aq, it is straightforward 
to integrate them out by diagonahzing the interaction terms. The resultant baryon 
interaction potential V is determined from f dt V = — ^on-shcu as 

y _ 277riV. 1 {{tm)' + pI + pI) (4(rM)^(p^ + pI) - uplpl) 

XMkk pIpI mrMYY - bu{rMnpl+pl) - u {p\+pl - {u + 2)plpl) ' ^ 

where u = (tr [f/(i)t;7(2)] )2 _ 

In addition to this potential V, we have another baryon interaction potential 
which comes from the mass term of in the matrix model action It is easily 

evaluated as 

■ iMl^triXr = ^M|k {{ny/2 + 2{Y,f) . (4.17) 

Here (^4)^ term is the mass term in the single-baryon Hamiltonian, so it does not 
contribute to the two-baryon interaction. On the other hand, the off-diagonal ele- 
ment is intrinsically the interaction between the baryons. The definition ( [4.7|) is 
computed as 

= -^r, tr[2r^(f/«)t[/(2)] ^4_^g) 
where i = 1,2,3, so we can write the potential energy explicitly as 



2 , pIpI ' ' ^ -^^2' 



1^4) +7-:rT^(^itr 

\' M) 



1627r 

Therefore, in total, the two-baryon interaction Hamiltonian is given by 
277riV, 1 ((rM)' + P? + pi)(4(rM)'(p? + pi)-«pfpi) 



(4.19) 



V 



AMkk p\pI \^{{tmYY - ^n{TM)\p\ + pD - ^ (pf + P^ - + 2)p?pi) 



, aat^mIk p\pI . ^ , 

+ 1627r [Tl,f^'^^' 



where u = (tr [[/'^^^'''t/*^^-'])^ — 4. The last term is the subtraction of the single-baryon 
Hamiltonians (the first term in ( p.l8| ), while the second term (mass term) in ( p.l8| ) 
cancelled already). 

Next, we shall evaluate this potential for a given quantum state of the two 
baryons, and show that there is a universal repulsive core for any baryon state. 

4.3 Universal repulsive core 

First notice that in the classical limit, i.e. Nc — )■ 00, the expectation value of any 
function (/(pi,P2)) with any given quantum state of the baryons approaches the 
classical value /(pi = p, p2 = p) where the classical p is given by ( |3.19| ). The 
equivalence between iVc —)■ 00 limit and h ^ limit is due to the fact that our 
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matrix action ( p.l| ) has an overall factor N^. Since any deviation from the classical 
value is 0{1/N^), while we are keeping only leading terms in the large Nc expansion, 
we in effect just need to put pi = p2(= p) in our potential ( [4.2UD .^^ We find 



V 



21tiN, 



(tr [f/Wt[/(2)])- 



4AMkk {tmY + 2p2 - 1 (tr [f/(i)tf/(2)])' p 



1627r 



(4.21) 



M) 



This is positive semi-definite. In fact, we can show that our potential ( [4.21| ) never 
vanish. In order to see ( |4.21| ) never vanish, note that it could vanish if and only if 
tr [C/(^^'^?7'^^)] = and tr [ir*f/(-'^)"l't/*^^^] = for any i. However the latter condition 
implies f/*^^)^f/^^^ oc 12x2? which contradicts with the former, therefore this is impos- 
sible. In this way, we see that ( [4.2 1| ) is positive definite and as a result, there is a 
universal repulsive potential (core) for any choice of two-baryon quantum states. 
If we expand ( [4.201 ) for (r^)^ ^ p^, we obtain a leading term 



V 



27'nNr 



64AMkk {tmY 

, xnmIk pIpI 



2 ^2 



M 



1627r 

For the large iVc, again we can put pi 

21tiN, , ^..M^....o^.^2 1 



6 (tr [t/Wtf/(2)])% [t/Wtf/(2)])- 



P| , Pi 

PI pi 



rj tr 



(4.22) 



P2 = p and obtain 



V 



4AMi 



KK 



(tr [[/«tf/(2)]y 



AiV,M3 



KK P 



16271 



r.2 )2 

Ml 



(r,tr [zr^'f/«t[/(2)])2 ^423) 



Again, we see the repulsive core. The universal repulsive potential scales as 1/r^ at 
r ^ p where r is the baryon separation. 

Let us consider the vacuum expectation value of the Hamiltonian ( [4.231) . The 
quantum state of the two baryons is specified by (7^, J^, np 



(1) ^» 



n 



with i = 1,2 which 

labels the two baryons. / (J) is the isospin (spin) of the baryon, while Up and nz are 
labels for excited baryon states ||16|. Explicit wave functions are given in [jl6|. For 
nucleons (|/| = \ J\ = 1/2), the spin/isospin wave functions are 



TT 



IJ 



bt) \pi) 



(4.24) 



^^This argument is valid as long as the classical value does not vanish and the wave function 
localizes for Nc ^ oo, which is in fact the present case for p and X'^. As for the spin/isospin 
encoded in U'--^^ and [/^^^ we cannot take the classical value as the wave function is not localized in 
the group space. 

^■^ Although the vector tm is in 4 spatial dimensions, once we take the VEV with the wave function 
of X^, it reduces to a 3-dimensional vector {i = 1, 2, 3). This is again because the leading term 
in the 1/Nc expansion is a classical value of X* which is zero. 
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for each nucleoli, so for our two-nucleon case the wave function is 

1 



71 



{r'U^%.j,{r^U^%,j, (4.25) 



where (Ji, Jiyh, J2) are the third components of the isospins and the spins for the 
nucleons, thus take values ±1/2. 

Now we obtain an explicit repulsive core for nucleons. At the leading order in 
l/Nci we can simply take the classical values for p and X^. As for the spin/isospin 
part, we can use the formulas given by (4.9), (4.11), and (4.18) of 



(tr[t/«tt/(2)])n = l+^-^IlPM, (4.26) 
(tr[^r^f/(i)tf/(2)] tr[«r^f/(i)tf/(2)]^^^^^^^^^^^^= + ]^jkjk ^J^^^^ ^J^^_^ij JkJk^^ 

These formulas are for a given isospin I\ and spin J\ for the nucleon labeled as 
A = 1,2. Finally, the vacuum expectation value of the potential ( [1.23| ) gives the 
central and the tensor forces, 

{V)i,,JM = Vb(r-) + SuVTif^ (4.27) 
with the standard definition 5*12 = 12J{fj Jg^j — 4J{ (with fj = ri/\r\), where 

This is the short- distance nuclear force obtained from our matrix model. We find 
there is a repulsive core of nucleons. The repulsive potential scales as 1/r^ for the 
inter-nucleon distance r, which is a property peculiar to the holographic model, as 
noted in [gg.i^ 

Finally let us compare this result with the soliton approach ||20| in which the two- 
baryon interaction potential \/soiiton^^ _ y^oiitonj^j^ _|_ ^^^y^ohton^^^ computed as, 

yc''-yr>^^[^,frjiJiJi)^J^, (4.30) 
^-o>-,.|/;/._^J,. (4.31) 



^^Thc nuclcon-nucleon potential normally consists of two different regimes. One is the long 
distance regime of order 1/Mkk and beyond ||2l|, where the four-dimensional exchanges of relatively 
light mesons, including pions of course, dominate nucleon interaction and induces attractive forces. 
The other is the repulsive core at short distances. One might naively think that the latter is not too 
relevant for formation of nuclei for example. However, nuclei are finely balanced systems with very 
small binding energy, and at least numerically seem sensitive to the details of the short distance 
repulsive core. 
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Compared to our ( |4.28| ) and ( |4.29| ), we can see that the structure of the core is the 



same.^^ Furthermore, the numerical coefficients are only different by factor 5/4 for 
the term proportional to im^iJi Vc°''*°"(r) and for the term I/r^"''*""^ (r ) . One can 
see easily that the difference in the last coefficient by factor 5/4 is due to the relation 
( p.25|) . More generally, recall that in the large one-boson exchange picture of 
nucleon-nucleon potential |^ , the latter two isospin-dependent structures are known 



to arise only from the second class of couplings, g^^^\ of section p.5| . It is the coupling 
squared that enter the potential, which explains the factor 5/4 increase relative to 
the soliton model. The isospin-independent central term arises, again in the large 
Nc, from iso-singlet vector exchanges, couplings for which belong to the ffist class, 
g'^^\ which explains the agreement. 

It is encouraging that the soliton approach and our matrix approach give qual- 
itatively the same, and also quantitatively similar result here as well. The large 
Nc computation is relatively stable in interpolating short distance and long distance 
regime. At the same time, numerically the increase is hardly negligible and may yet 
prove to be a big difference in the end, when we have a good control over subleading 
1/Nc and 1/A corrections in the model. Whether or not this will favor the current 
D4-D8 model in simulating QCD is unclear for now, however, since some of known 
subleading l/Nc corrections appear at comparable magnitude ii Nc = 3 is used and 
it is not known how to catalog all such corrections. 

Before closing, we would like to make comments on some of existing computations 
of nucleon-nucleon potential. First of all, we again emphasize that we only computed 
the repulsive core at very short distance. As we mentioned above, this short-distance 
repulsive core was computed by the solitonic methods in Ref. |^ and also in Ref. [P^] , 
which differ from our result only by a numerical factor. The fact that there is a 
universal 1/r^ behavior is a clear and unambiguous prediction of holography and 
needs to be compared against other less speculative approaches, possibly via lattice 
QCD simulations or by experiments. 

Of more practical interest to most of nuclear community would be the long- 
distance behavior, which should include an attractive iso-singlet channel. From the 
holographic viewpoint of D4-D8 model, this long-distance behavior was found in 
Ref. where the authors computed one-boson exchange potential, with coupling 
constants computed precisely using the sohtonic picture of the baryon, and indeed 
found the long-distance attractive potential with right qualitative behaviors, such as 



allowing deuteron wavefunction numerically p5[. Our matrix model, in its current 



form, cannot address this long distance behavior, unfortunately. 



"'^^It is also interesting to note that, at the stage of the Hamihonian, the second term of (4.20) 
completely coincides with what is called ijp'^^^)) in |^ in its structure. This part basically gives 
the tensor force. The difference in coefficients in the tensor forces are due to the relation ( ^.25| ). 
Nevertheless, since /f^'^'^') came in from a generalized Osborn's formula, it is quite interesting 
that our matrix model can reproduce quite easily the formula. 
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There exist extensive literatures on the matter of nucleon-nucleon potential, 
which are too numerous and too diverse to list. See for instance Ref. for com- 
putations based on the chiral perturbation theory. In the face of such extensive 
previous work, one may ask why we should again try to compute nucleon-nucleon 
potential. The first reason is that at least the top-down models such as ours have 
very few adjustable parameters and should be very predictive and precise. This is 
not to say, of course, that the latter would be more accurate. On the contrary, such 
top-down holographic computations are less likely to be near the correct low energy 
QCD physics. 

What we aim to build is not the most phenomenologically correct model, but 
rather those which can give us insights and results from top-down approach and that 
cannot be given by any of conventional phenomenological methods. The universal 
repulsive core above is one fine such example. The conventional chiral pertur- 
bation, for instance, generates type repulsive core from exchange rho mesons 



via tensor couplings as well as 1/r from omega meson exchange whose precise 
behavior, however, should be taken with a grain of salt at short distance below the 
mass of the nucleon. In practice, one further introduces other short-distance cut-offs 
to better simulate nature. Whether or not the holographic repulsive core found 
by above holographic computations represents a better solution to this short-distance 
treatment of nuclear force and a more faithful image of real QCD is a matter to be 
settled eventually by lattice QCD simulations or by experiments. 

5. Conclusion 

We have derived a U{k) matrix model which describes fc-baryon systems at 

short distance, by considering open string theory on the wrapped baryon vertex D- 
branes embedded in the D4-D8 model of large Nc holographic QCD. With this matrix 
model, we computed the holographic size and the spectrum of the baryon {k = 1), 
and also short-distance nuclear force {k = 2). The latter exhibits a repulsive core, 
which has been quite important in nuclear and hadron physics for long years. 

This model complements the Yang-Mills soliton picture of the holographic baryon 



Ig, |16|, in that it is capable of addressing short-distance behavior such as the 
crucial hadronic size estimate of the baryon. Recall that the size found in these 
studies turned out to be comparable to the (appropriately warped) string scale, 
which is too small to be justified by the Yang-Mills method used there. Our new 
matrix model is trustworthy well below the string scale, in the opposite end of length 
scale, and represents an opportunity to check whether the soliton picture gave us the 
correct estimates or not. As we summarize shortly, the matrix model estimate gave us 
essentially the same size up to a numerical factor of (5/4)^/'^ ~ 1, showing that stringy 
correction does very little to correct the soliton picture. This is very fortunate for 
the underlying D4-D8 holographic QCD since, at least for small number of baryons. 
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the soliton picture has generated many predictions that agreed with nuclear data 
[p!7| , p!8| , |19|, This is by far the best evidence we know of that justifies study of 
baryons in holographic QCD, despite the latter's huge mass and small size in the 
large Nc and large A limit. 

A distinguished character of our matrix model is its simplicity. In particular, the 
form of the matrix model does not depend on the number of the baryons described 
— simply changing the rank of the matrix allows one to treat larger number of 
baryons. By generalizing the method we did in this paper, we can handle general 
A;-body systems. This is significant simplicity compared with the soliton picture, 
where handling general fc-solitons are quite painful due to the significant increase of 
the number of moduli parameters. 

Therefore it is natural to apply our matrix model analysis to more general k- 
body systems. In nuclear physics, the role of 3-body forces is very important. For 
example, for few-body nuclear bound states (hght nuclei), the three-body forces are 
known to play a crucial role. Three-body forces in the soliton picture in the D4-D8 
system were derived in ||3^. However, in that paper the isospins and spins were 
treated classically due to the complications induced by the increase of the number of 
moduli parameters, so the analysis was not complete. Applying our matrix model, 
it is possible to treat generic spins and isospins quantum- mechanically ||3^. It would 



be interesting to see how the quantum spin and isospin effects modify the previous 
soliton picture results p7 . 



It would also be interesting to consider a large k limit, where the states are 
similar to extremely heavy nuclei or a core of neutron stars. The reader may wonder 
the effect of a back-reaction in that limit. Taking the large k limit certainly cause the 
back-reaction of the multiple D4'-branes in D8-branes in D4-D8 system. However as 
we have seen, in our system, we have no supersymmetric leftover, and as a result, 
there is a repulsive forces between D4'-branes. So the D4'-branes cannot approach 
closer than p = 0{l/\/X) ~ I'f , which is too large compared with p ~ {I'fy for 
Maldacena's decoupling limit |^ in the If^ — ?■ limit. Therefore the back-reaction 
is expected not so strong enough to make a back-reacted gravitational throat. 

Going back to the results for k = 1, 2 in this work, we note that the results agree 
with the previous soliton approach qualitatively, and in some case quantitatively 
also. The central isospin-independent part of the repulsive core is one example. The 
holographic size squared of the baryon also turned out to be larger only by a fraction 
which, considering the vast separation of scales between the two approaches by a 
factor of A^/^, seems pretty innocuous. This tells us that physics of baryon in D4-D8 
holographic QCD is relatively stable against corrections. This stability of the baryon 
physics against short distance corrections was previously anticipated and argued for 



^^This argument could break down for finite if^ or in supersymmetric cases and we may consider 
a holographic dual of the multi-baryon system in the fc — ?> oo limit. It leads to "holographic nuclei" 
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on the basis of an approximate supersymmetry and is borne out in our new 
matrix model approach. Without such stabihty, study of baryons in holographic 
QCD would have been very difficult and cumbersome. 

On the other hand, the new results also show nontrivial change of key observables, 
which cannot be ignored. From the viewpoint of low energy effective theory in terms 
of mesons and baryons, the change can be summarized in the large limit as 
a universal and multiplicative increase by (5/4)^/^ of certain meson- nucleon-nucleon 
couplings. These are couplings denoted as g^^^^ in section 3.5, and include the leading 
axial couplings to pions, all isotriplet axial vector couplings of the operator dimension 
four, and the dimension-five derivative couplings to isotriplet vector mesons. In 
particular, the change came about not as subleading l/Nc corrections, but a change 
in coefficients of the leading terms. 

The reader may wonder whether the matrix model can describe only the short 
distance physics, such as the repulsive core of the nucleon. As we argued in section 
3.5 and reiterated above, the model is capable of addressing the long-distance physics 
of soft-meson processes, via the single quantity p^, so in that indirect sense also tells 
us much about baryons far away from one another. An entirely different question is 
whether we can use the matrix model directly to compute interaction, as we did in 
section 4, between well-separated baryons without resorting to intermediate mesons. 
This is particularly important when we wish to consider large k physics where meson 
exchange viewpoint would become quickly intractable. 

The above matrix model does take into account the dual warped geometry of 
the pure QCD background, so it already knows about effect of glues at energy scale 
Mkk- Yet, how to emulate the long-distance physics due to exchange of pions and 
other light mesons directly by manipulating the current matrix model is a more 
challenging problem. For example, in a similar problem of DO-brane interactions, 
it was the one-loop effect in the open string side that emulated the long distance 
graviton exchange between them. This agreement was explained by the combination 
of supersymmetry and the worldsheet channel duality [^. In the present case, we 
must emulate an exchange of D8-D8 open strings by a matrix model of D4'-D4' and 
D4'-D8 open strings. This appears to be a qualitative different problem from the 
case of DO interactions. Absence of unbroken supersymmetry probably make matters 
a bit worse also. 

Nevertheless, we believe that this can be achieved in the matrix model with 
judicious insertions of effective operators and by considering quantum effects. Once 
this is done, on the other hand, our matrix model may prove to be much more 
versatile than the previous instanton soliton picture of baryons, since it can be more 
easily generalized to many nucleon systems. Nuclei with generic atomic number k and 
perhaps other dense matter systems would be more accessible that way. Perhaps we 
can eventually construct a U{k) matrix model, where attractive long-distance forces 
and repulsive short-distance forces are carefully balanced and produces,say, helium 
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nucleus and carbon nucleus. 

In this context, we finally comment that we only made use of bosonic fields in 
this work, even though there should be fermionic partners in D4'-D8 matrix theory. 
For this work, we did not need loop computation, which justifies this truncation. 
For a more complete matrix model capable of directly dealing with baryons at large 
separation, we probably need to consider fermionic fields as well, which is beyond 
the scope of this paper. 

Eventually we hope to be able to compute other interesting quantities in our 
matrix model, such as currents associated with chiral symmetry, energy-momentum 
tensor and equations of states from that. In particular, equations of states with 
large baryon number k are important to understand the dynamics of astrophysical 
compact starts, such as neutron stars. 
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